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Abstract. We introduce multi-scale filtering by the way of cer- 
tain double convolution systems. We prove stability theorems for 
these systems and make connections with function theory in the 
poly-disc. Finally, we compare the framework developed here with 
the white noise space framework, within which a similar class of 
double convolution systems has been defined earlier. 
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1. Introduction 

A wide class of causal discrete time-invariant linear systems can be 
given in terms of convolution in the form 



'1-1) Z/n = hn-mUm, ^ = 0,1, 



m=0 



where is the impulse response and where the input sequence (um) 
and output sequence belong to some sequences spaces. The Z 
transform of the sequence (hn) 

oo 
n=0 

is called the transfer function of the system, and there are deep rela- 
tionships between properties of h and of the system. For instance the 
system will be dissipative in the sense that for all £2 inputs it holds 
that 



2 

n| ) 



E £ E 

n=0 n=0 

if and only if h is analytic and contractive in the open unit disc D. The 
function h is then called a Schur function. Similarly, the system will 
be ii — I2 bounded if for every l\ entry, the output is in I2 and there 
is a M > independent of the input such that 

\ 1/2 

vn=0 / n=0 

As is well known, the system is i\ — £2 bounded if and only if h be- 
longs to the Hardy space H2(D). Systems of the form (11. ip . the Hardy 
space 112(115) and Schur functions have been generalized to a number 
of situations in the theory of A^-dimensional systems and beyond; see 
for instance the works [T], [TU] [TT], [T^, and the references therein. 

Another generalization of systems of the form (II. ip occurs when and 
Vt^ri ctre not complex numbers, but belong to some space with a product, 

say 



'1.2) yn = ^hn-m'^Um, n = 0, 1, 

m=0 



Of special interest is the case where ^ir is a convolution as, for instance in 
[1], and the system is then called a double convolution system. Therein, 
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the first named author together with David Levanony considered an ex- 
ample of such a double convolution system, when both and are 
random variables, which belong to the white noise space, or more gen- 
erally to the Kondratiev space. The product hn-mUm in (H-ip is then 
replaced by the Wick product. The Wick product takes the form of a 
convolution with respect to an appropriate basis, and we have an ex- 
ample of a double convolution system. Using the Hermite transform, 
one can define a generalized transfer function, which is a function ana- 
lytic in z and in a countable number of other variables (these variables 
take into account the randomness). The white noise space setting is 
reviewed in the last section of this paper, with purpose the comparison 
between the present paper and 

In the present work we study another type of double convolution sys- 
tem, which arises in the theory of multi-scale systems. We use the 
approach of the second named author presented at the Mathematical 
Theory of Networks and Systems conference in 2006 in Kyoto, see |22|, 
to define the multi-scale version of the systems fll.ll) . Let 

and let = ^ 

be the corresponding automorphism of D. Following [22j, consider the 
map 

cz + a 

where / is analytic in D. Then T^pf is also analytic in the open unit 
disc. Let 

oo oo 

/W = E/"^" and {T^fKz) = Y,U^^z^ 

n=0 n=0 

be the Taylor expansions of / and Ti^f respectively. Let 
N={1,2,...} and No = {0, 1, 2, . . .} . 

The map which associates to the sequence (/«)neNo the sequence (/«,(^)neNo 
is called the scahng operation (see the precise definition in Section [2]). 
Consider now a subgroup F of 5*^7(1, 1), which represents the scales we 
will use to study the signals and systems. One associates to the se- 
quence (/n)„gNo its scale transform (/n(7))neNo,7Gr, which is a function 
of n G No and 7 G F and where we have written fn{j) for fn,-y In the 
case which we will consider, F will be indexed by Z^, but the resulting 
setting is quite different from classical ND theory. 
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The scale transform is the starting point of our approach (initiated in 
[22]) to multi-scale analysis in discrete time. In opposition to wavelets, 
we propose a transform which has on the same level both the time 
and the scale aspect. Let us now elaborate on the differences of our 
approach to wavelets. Recall that in continuous time the wavelet trans- 
form of a signal / is defined by 



Wf{u, 



where s is the scale parameter and ip is the mother wavelet; see 
pp. 78-79]. The discrete time version of the transform is obtained, for 
discrete scales s = , by discretizing the above integral as in 



Wf[n,a^] = Y,fkt 



k-n,j> 

where ipnj = '^4' (^) cind fk = f{k), assuming a sampling period 
normalized by 1 (see [211 pp. 88-89]). The decomposition provides 
an appropriate mathematical tool for signal analysis. In particular, it 
makes it possible to extract the components of a given discrete-time 
signal at a given scale on a discrete grid. However, the question of 
defining the scale shift operator (dilation) for purely discrete-time sig- 
nals is dodged somehow. Another point of departure from our approach 
is that the wavelet transform (either continuous- or discrete-time) has 
only one convolution as compared to the double time and scale convo- 
lution considered in the present work (see equation (11. 3p below). 

We define linear systems as expressions of the form (11.21) : 

n 

Un ^ ^ ^71—771 ^ ^771; 

m=0 

where -k denotes the convolution in V. Thus, 

(1.3) ?/„(7) = ^ I ^/i„_„(7 0(^-i)u„((^) J , 7 e r. 

ni=o Viper / 
When r is trivial we recover (11. ip . 

There are parallel and analogies between the theory of linear stochastic 
systems presented in [1] and the theory developed here. These parallels 
are pointed out in the sequel, and serve as guide and motivation for 
some of the proofs in the present paper. Still, there are some differ- 
ences between the statements and the proofs of the stability theorems 
in [1] and the proofs given here. These differences are pointed out in 
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the text. We note that a general theory of double convolution systems 
is in preparation, [S]. We also note that some of the results presented 
here have been announced in [7|. 

We now turn to the outline of this paper. It consists of 10 sections 
besides the introduction. In Section [2] we review the approach to dis- 
crete multi-scale systems presented in [22]. In Section [3] we define the 
systems which we will study in the paper, and the related notion of 
transfer function. Section H] is of a review nature. We discuss the 
trigonometric moment problem and related reproducing kernel Hilbert 
spaces of the type introduced by de Branges and Rovnyak. In Section 
[5] we consider the case where the sub-group has one generator and is 
infinite. We use the classical one dimensional moment problem to asso- 
ciate to the Haar measure of the dual group a uniquely defined measure 
on the unit circle. This allows us to use function theory on the disc 
and on the bi-disc. In Section [6] we review some deep results of Mihai 
Putinar, see [23], on positive polynomials on compact semi-algebraic 
sets and their use to solve the trigonometric moment problem on the 
poly-disc. In Section [7] we consider the case where F is no more cyclic 
but has a finite number of generators. Although the results in Section 
[5] are particular cases of the ones in Section [7] we have chosen to present 
both for ease of exposition. The next three sections consider stability 
results: BIBO stability is considered in Section [8], dissipative systems 
are studied in Section [9] and Section [10] is devoted to ii — £2 stability. 
In the last section we review the white noise space setting, and present 
a table and some remarks, which point out the analogies between the 
setting in [1] and the present work. 

Acknowledgements: It is a pleasure to thank Professor Mihai Puti- 
nar for explaining to us the solution of the moment problem in the case 
of the poly-disc. 

2. Scaling operator for discrete- time signals 

We briefly summarize the approach to multi-scale systems presented 
in • We first note the following: If F{s), ^ 0, denotes the 

Laplace transform of a continuous-time signal f{t), t ^ 0, then, for any 
a = 1/P > 0, y/aF{as) is the Laplace transform of f{Pt). Therefore, 
time scaling has a similar form in the frequency domain. As opposed 
to the continuous-time case, time scaling is not clearly defined in the 
discrete-time setting. Nevertheless, the preceding remark is the key 
step to define a scahng operator for discrete-time signal. Consider the 
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Mobius transformation 

which maps conformally the open right half-plane C+ onto the open 
unit disc. To recall our definition of the scaling operator (see [22j and 
also [8]), we note that the scale shift in C+, 

Sa'- s ^ Sa{s) = as, a > 

translates in the unit disc, via Gels), into the hyperbolic transformation 

(e^^ + ae~'^) z + (1 - a) 



(2.1) ^{a}iz) = {GeoSaoGe')iz) 



'1 - a)z + (e-*^ + ae*^) 



Any such transformation maps the open unit disc (resp. the unit circle) 
into itself. Now, the most general linear transformation which maps 
the open unit disc (resp. the unit circle) into itself has the form 

2.2 72; = , 7i — 72 = 1. 

If |3?(7i)| > 1, then the transformation is hyperbolic [15] and it takes 
the form 

(2.3) = -, > 0, 

7(2;) - ^2 

where = V [^(71 1+^^(71 ) a _a; 

7j -y* -y* 

points. The constant is called the multiplier of the transformation, 
see [IHl p. 15], and is given by 



3?(7i) - v/[3^(7i)F-l 



3fJ(7i) + vmiiw - 1 

Noting that |^i| = |^2| = 1, one may rearrange (12. 3p to obtain 

- XZe'^-'jiz) X., - X*e'^-'z 

= a 

1 + 6*^^7(2;) 1 + e^^^'z ' 

where e*^^ = Dividing both sides of this equality by |A^| and setting 
e*^^ = we recover (12.11) up to a rotation: 

(2.4) e^«^7(^) = {Ge, o S^^ o G,;)(e^«^^). 

Any hyperbolic transformation 7 of the form (12. 2p is thus (conformally) 
equivalent to a scale shift 5*0, in C+. 
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In the sequel, we will be interested in Abelian subgroups of hyperbolic 
transformations, but the remainder of this section deals with general 
linear transformations 

(p(z) = — with ad — be — 1, 

from the open unit disk onto itself. To each such transformation, we 
associate (in bold letters) 



and we define 



cz + a 



Lemma 2.1. Let cpi and ^2 belong to SU{1, 1). Then 

(2.5) T^,oT^^=T^^^^, 

and in particular for every (p G SU{1, 1) it holds that 

(2.6) T^-. = {T^)-\ 
Proof: We have 

1 1 

C2Z + d2 Ci(p2{z) + di 

But 

11 11 



C2Z + d2 ciip2{z) + di C2Z + d2 ci^j^ + di 



C2Z+d2 
1 



(cia2 + diC2)z + C162 + did2 ' 
which ends the proof since the second row of 

/ ai bi\ f 02 62 

\Cl dij \C2 d2 

is equal to 

(cia2 + diC2 C162 + did-^ . 



□ 



Let / be analytic in the open unit disc. Then T^^f is also analytic 
in the open unit disc. The mapping induces a mapping from the 
space of sequences coefficients of power series of functions analytic in 
a neighborhood of the origin into itself: if f{z) — Yl^=o z^^n is the 
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power series expansion at the origin of the function / analytic in D, 
then 

oo 
n=0 

In view of f l2.6p we have: 

Proposition 2.2. Let e SU{1, 1), and let (a„)„gNo be a sequence of 

complex numbers such that hmsup„ \an+i\^^^^~^^^ < 1- Then there 

exists a sequence of complex numbers (&„)„gNo such that 

hmsup|6„+i|i/("+^) < 1, 

n >oo 

and 

a^i bn^(p-, ri 0, 1, . . . 

Proof: Let a{z) = Yl'^=o^nZ'"- It suffices to define a series 6„_(^ via the 
formula 

oo 

^&n,<^2;" = (T<^-ia)(z), 

n=0 

and use ^M)- □ 



Theorem 2.3. The operator Tip is unitary from H2(D) onto itself with 
norm equal to 1. It is also continuous from Hoo(ID)) into itself with norm 
equal to l/{\d\ — |c|). 

Proof: First recall the formula 

1 — zw* {cz + d){cw + d)* 

where z, w are in the domain of definition of (p. Furthermore, recall 
that a function / defined in D is analytic there and belongs to H2(]D'), 
with ||/||h2(b) < 1, if and only if the kernel 

(2.8) -^-f{z)f{wr 

1 — zw 

is positive in D; see for instance [U Theorem 2.6.6]. We now compute 
for 

Mz,w) = - {T^f){z){T^f{w)r 

1 — zw* 
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for z,w G D. Using (12. 7p we can write: 

1 



(p{z)(p{w)*){cz + d){cw + d)* 

^ f{^{z))mw)y 



The kernel 



{cz + d){cw + d)* 
1 



{cz + d){cw + d)* 
1 

1 — (p{z)ip{w)* 

1 



f{^{z))f{^{w)r 



f{v{z)f{v{w)r 



{l-v{z)v{wY] 

is positive in D since the kernel fl2.8p is positive there. It follows that 
A{z,w) is positive in the open unit disc, and thus, by [1, Theorem 
2.6.6], the function T^{f) belongs to H2(D) and has norm less or equal 
to 1. To prove that the norm is indeed equal to 1 we use (12.51) and 
(12.61) . which imply that 

l<||T<,||-||T<,-i||, 

which, together with the fact that both T^p and Ti^-i have norm less 
that 1 implies that \\T^p\\ = 1. We now show that is unitary. Let 
/ G H2(ro). We have: 

H2(D) = \\Ti^-^T^{f )\\H2m 

< \\T^if)\\u2io) 

< II/I|h2(d), 



since both and T<^-i are contractive. It follows that is unitary. 

The second claim is easily verified. □ 

We therefore associate to the signal x = {xn}nmo ^ ^2 the signal 
{xn <^}neNo indexed by No x SU (1, 1). In the sequel, we will simplify 

(p&SU{l,l) 

the notation by writing Xn{^p) in place of x„^<^. For fixed m, {xm{^p)}, 
with (p G SU{1,1), represents a scale signal, that is, the observation 
of the signal at time m through the scales cp G SU{1, 1). In the 
sequel, we will consider (by convention) the zooming as corresponding 
to the "positive" scales. 
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Definition 2.4. The scale-causal projection of {xm{f)}, ^ G SU{1, 1) 
is given by the restriction of {xmi'^)} to the scales ip for which the 
multiplier is strictly less than one: a^p < 1. 

Definition 2.5. Given a discrete subgroup T of SU{1,1) we denote 
by the set of transformations consisting of the identity and of the 
scales if for which the multiplier is strictly less than one: < 1. The 
system (11.21) will be scale-causal if the elements hn G £2(r+). 

Given 7 and ip two elements of F, we will say that 7 succeeds ip and 
will note <y9 ^ 7, if 7 o g r_|_ that is: 

V9 ^ 7 <^==^ a^o^-i ^ 1. 

Proposition 2.6. The relation ^ defines a total order in T. 

Proof: Since we assume that F is Abelian, all the transformations 
must have the same fixed points. The parameters and 9^ in (12.41) 
are therefore constant. The proof then follows upon noting that the 
multipher a^o^ is given by: □ 

With this order we obtain a bijection 

7 ^ ^(7) 

between F and Z, and one can identify ^2(r) and i2{'^) and £2(r+) and 
^2 (No). 

Remark 2.7. Using the isomorphism we introduce the following defini- 
tion: 

Definition 2.8. The function ^(7) from F+ into C has finite support 

N(u) = max {^(7) such that u{'y) 7^ 0} < 00. 
The support of the function u is the interval [0, N{'u)] C Nq. 

The results of this section remain valid if we replace the Hardy space 
H2(]D)) by its vector-valued version H2(]ri')®7i, where Ti is some Hilbert 
space. One can define in particular the scaling of random sequences 
when ?i is a probability space. 

Remark 2.9. In p2t equation (17)] another kind of systems are consid- 
ered, with only one convolution. The main issue there is the notion of 
scale-invariance in a stronger form, which will not be considered here; 
see also [6] for related work. 
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Finally, we note that one could consider systems non-causal with re- 
spect to n, that is of the form 

Thus, there are really four possibilities for the various stability theo- 
rems we present, depending on whether we have time causality or not, 
and scale-causality or not. In this paper we only give part of all possi- 
ble results. 



3. Discrete-scale invariant systems and signals 

The scaling operators T^^ form a group of operators from the Hardy 
space H2(]D') onto itself. From now on, we discretize the scale axis and 
restrict <^ to a discrete subgroup F of SU (1,1). We will take F Abelian 
(cyclic) and consisting of hyperbolic transformations, and we denote 
by F its dual group. Recall that F is formed by the set of functions 

0" : F — i> T such that a^i) = 1 and V 7, ip, (7(7 o ip) = a{'y)a{ip), 

where l stands for the identity transformation. The elements of F are 
called characters of the group F (see [16]). We denote by fl the Haar 
measure of F, which is compact by the Pontryagin duality We 
recall the definition of the Fourier transform on F and of its inverse: 

7er 

x(7) = j^x{a)a{-f)d'i2{a). 
The Haar measure djl is normalized so that Plancherel's theorem holds: 

ll/lll(r) = E I/WI' = i \f(^)\"dl^{cr) = \\f\\li,ny 
See [H Theorem 8.4.2 p. 123]. 

Definition 3.1. A signal will be a sequence {un{-)}nef>io of elements of 
£2(F), and such that the condition 

(3.1) sup l|M„(-)ll^2(r) < 00 

n=0,l,... 

holds. 

A scale-causal signal will be a sequence {Mn(')}neNo (^f elements of 
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£2(r+), and such that the condition 

(3.2) sup ||un(-)ll^2(r+) < oo 

n=0,l,... 

holds. 

In the sequel, we will impose the following stronger norm constrains on 
a signal, besides (13.11) or (13.21) . namely: 

oo 

(3.3) ^ ||M„(-)|||(r) < oo, 

n=0 

or 

oo 

(3-4) ^||M„(-)||£2(r) < oo, 

n=0 

and similarly for scale-causal signals. 

We note the following: a dissipative filter cannot be effective at all 
scales. At some stage, details cannot be seen. These intuitive facts are 
made more precise in the following proposition. 

Proposition 3.2. 

(1) Assume that the supports of the Un are uniformly hounded. Then, 
(13.31) is in force. 

(2) Assume that the support of Un is infinite for all n. Then, the sum 
on the left side of ( 13. 3p diverges. 

Proof: Let be such that the support of all the functions 7 ^ Uni'y) 
is inside [0,A^]. Then, 

00 00 N 

n=0 n=0 £i(7)=0 

N 00 

e{-y)=0 n=0 
< ^lkn||^2(No); 

since (see Theorem 12. 3 p the maps T^ are unitary from H2(]I3) onto itself. 

The second claim is proved similarly. □ 

An example of (m„) satisfying Condition (1) of the preceding propo- 
sition has been presented in the paper [22], where the corresponding 
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group r is Fuchsian. This was used therein, to define the scale unit- 
pulse signal. A similar condition was also considered by P. Yuditskii 
[25] in the description of the direct integral of spaces of character- 
automorphic functions. 

Definition 3.3. An impulse response (resp. a scale-causal impulse 
response) will he a sequence {/in(')}neNo of elements of i2(X) (resp. of 
i2{^+)) such that for every n G Nq, the multiplication operator 

(3.5) -^h„ : u^hn-ku, 72 = 0,1,... 

is hounded from £2(r) into itself (resp. from £2(1^-1-) into itself) and 
such that 

(3.6) sup ||/i„||£2(r) < 00 {resp. sup ||/i„||£2(r+) < oo)- 

n=0,l,... n=0,l,... 



The systems that we consider here are defined by (11.21) . that is, by the 
double convolution (II. 3p . that we recall below: 

(3.7) yn{l) = I X] ^"-'"('y ° ^~^)Xm{.8) j . 

■mez \5er / 

In view of (13. 6p the series 

00 

(3.8) i7(z,a) = ^^"/i„(a) 

n=0 

converges in the L2((i/i) norm for every 2; G D. Taking the Fourier 
transform (with respect to F) of both sides of (11.31) we obtain 

where the equality is in the L2((i/2) sense. Taking now the Z transform 
we get 

(3.9) Y{z,a) = H{z,a)U{z,a), 
where 

00 00 

(3.10) Y{z,a) = Y,^"yn{cr) and ?7(z, a) = z"S„(a), 

n=0 n=0 

and where, for every 2; G D the equality in (13.91) is /i-a.e. 

The function H{z, a) can be seen as the transfer function of the discrete- 
time scale-invariant system. Formula (13. 8p suggests to define and study 
hierarchies of transfer functions, for which the functions hn depend on 
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cr in some pre-assigned way (for instance, when they are polynomials 
in a), or when the function H{z,a) is a rational function of z or of 
a. In the next two sections, under the hypothesis that the sub-group 
r has a finite number, say p, of generators, we will associate to the 
system (11.21) an analytic function of p + 1 variables, which we will call 
the generalized transfer function of the system. 

4. The trigonometric moment problem 

We first gather some well known facts on the trigonometric moment 
problem in form of a theorem. 

Theorem 4.1. Given an infinite sequence . . . , to, ti, . . . of complex 
numbers such that 

t^n = t*n, n = 0, 1,..., 
there exists a positive measure dv on [0, 2tt) such that 

tn= [ e-''''du{9), neNo, 
Jo 

if and only if all the Toeplitz matrices 

"^N (j^n—m)n,m=0,...,N 

are non-negative. 



See for instance |20l Theorem 2.7 p. 66]. The measure is then unique 
(when normalized). We also recall that the sequence (t„) and the mea- 
sure du are related by 



to + 2j2t^z-= / -^Me) 

„_i Jo e — z 



n=l 

and thus the function 



oo „ 

$(^)=to + 2Vt.z"= / 

n=l 



e'' + z 
e*^ - z 



du{d) 



is analytic and has a positive real part in the open unit disc. Using 
Stieltjes inverse formula, one can recover v from $ via the formula 

lim / Re{^{re'^)}dd = u{h) - i^(a_). 



where we assume that v is right continuous. We also recall that the 
function 



2{1 -zw*) Jo {e^^ - z){e''^ - w)* 
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is positive for w G C \ T. We denote by the associated repro- 

ducing kernel Hilbert space when z and w are restricted to the open 
unit disc. The following result has first been proved by de Branges and 
Shulman; see pLSj. In the statement, H2(c?z/) denotes the closed linear 
span in "h^idv) of the functions for m > 0. 

Theorem 4.2. The space consists of the functions of the form 

h{z) = / ^.^ ^ du{t), h G H2{du), 
Jo e** - z 

with norm 

\\h\\c+i'!>) = \\h\\u2{du)- 
We now recall some results on the structure of the space H2((iz/). 

Theorem 4.3. Assume that Ji2{di>) ^ L2((iz/). Then H2((ii/) is a 

reproducing kernel Hilbert space, and its reproducing kernel is of the 
form 

A{z)A{w)* - B{z)B{w)* 
1 — zw* ' 
where A{z) and B{z) are functions analytic off the unit circle. 

Proof: We assume that H2((iz/) ^ L2((iz/), and let Hq G L2((iz/) 
H2((ii/). Let « G C \ T be such that 

ho{e)du{d) ^ ^ 

e*^ — a 

Let p be a polynomial; then R^ip G H2(c?z/), where 

(fi„P)(.) = P<£)^. 

z — a 

Then 

and therefore we obtain 



p{a) 



e^" — a 



2n 



ho{e)du{e) 



a 



Therefore the map p ^ is continuous on the polynomials, and 
extends to a continuous map to H.2{di>). Therefore H2((iz/) is a repro- 
ducing kernel Hilbert subspace of L2(c?z/). The proof is then finished 
by using ^ Theorem 3.1 p. 600]. □ 
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5. The case of one generator 

In this section, we consider the case of a cychc group F, generated by 
a hyperbohc transformation 7q G SU{1, 1). Any transformation in F 

is thus of the form 7™ = 70 o . . . o 70, m G Z. 

m times 

Theorem 5.1. There exists a positive measure dviO) on [0, 27r) such 
that 

(5.1) lG{^'^)d4i{a) = e'^^^d^ie), m G Z. 

Jt Jo 

Proof: We use Theorem 14.11 Let 

tm = J^cri^oTdJlia), m G Z. 

Since |ct(7o)| = 1 for all cr G F we have 

/^(7o)V(7o"^)rf/2(a) = (a(7-),a(7^))L,(.^), 
and therefore all the Toeplitz matrices 

= {tl-m)i,m=0,...N 

are non-negative. It follows from Theorem 14.11 that there exists a 
uniquely defined measure dv such that 

tm= ! e-'''^'dv{e), m = 0,1,2,..., 
Jo 

and hence we obtain (15.10 . 

□ 



Remark 5.2. The proof of the previous theorem formalizes the intuitive 
idea that one can make the "change of variable" 

Theorem 5.3. The linear map I which to o"(7^) associates the func- 
tion : 

(5.2) I(a(7o"^)) = mGZ, 
is an isomorphism from Li2{d'fl) into h2{di'). 
Proof: For a function / of the form 

M 

(5.3) /((t) = ^c„(t(7^) where N,MeNo and c„ G C, 

-N 
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we have 



n,m=-N,...,M 

2tv 



n,m=-N,...,M "^^ 

/ I 5^ Cne''''\''du{e) 

Jo AT 



-N 
\'L2(du)- 

The result follows by continuity since such / are dense in L2 (ci/i). To 
verify this last claim we note the following: By Plancherel's theorem, 
the map from ^2(r) onto h2{d'j2) which to the sequence which consists 
only of zeros, except the ?7,-th element which is equal to 1, associates 
the function (7(70)'^, extends to a unitary map. □ 

We will be interested in particular in the positive powers of 7g, which 
correspond to zooming (we consider that the multiplier of 'Jq, i.e. the 
associated scale Q^'-yQ 5 IS less than 1). 

Definition 5.4. We denote by Ji2{djl) the closure in Li2{d'fl) of the 
functions (7(70)", n = 0,1,2,.... Similarly, we denote by ii2{di') the 
closure in L2((iz/) of the functions z"' , n = 0,1,2, .. .. 

Note that it may happen that h2{d'fl) = H2((i/x). 

Following [1] we introduce the next definition. 

Definition 5.5. The map I will be called the Hermite transform. 

Recall that T is compact and therefore 
(5.4) L2 (c//i) C Li((i/i). 

In general the product of two elements / and g in L2 (cZ/i) does not 
belong to L2 ((i/i), and one cannot define l{fg), let alone compare it 
with the product l{f)I{g). On the other hand, we will need in the 
sequel only the case where at least one of the elements in the product 
fg defines a bounded multiplication operator from L2 (c?/i) into itself; 
see Definition 13.31 and the proof of Theorem 19.31 for instance. This is 
exploited in the next theorem. 

Theorem 5.6. Let f G Ij2{d'fl) such that the operator of multiplication 
by f defines a bounded operator from L2 ((i/i) into itself. Then for every 



18 D. ALPAY AND M. MBOUP 

g in L2((i/i) it holds that: 

(5.5) lUg) = i{m9)- 

Proof: We note that the muhiphcative property (15.51) holds for / and g 
of the form (15.31) . To prove the theorem we first assume that g is of the 
form (15.31) . and consider a sequence (p„) of elements of the form (15. 3p . 
converging to / in the L2 (ci/i) norm. The function g is in particular 
bounded, and so fg G L2 (ci/i), and we can write: 

Wfa -Png\\i.2{dri) < K\\f - pnWuidii), 

where K > is such that \g\ < K. Thus, fpn tends in li2{d'fi)- norm 
to fg. 

The function I{g) is bounded, and I(/) G L2 (c?z/). Therefore: 

\\I{fg)-limg)\\uidu) < \\I{fg)~I{pn9)\\uidu) + 

+ ||I(p„c7)-I(/)I((7)||L,(d.) 

= Wfg -Png\\L2idp)+ 

+ \\i{pM9)-mm\\w 

< Wfg -Pn9\\L2{dn)+ 
+ Ki\\pn - gWLiidp), 

where iiTi > is such that |I(5')| < Ki, and where we have used that 

I{Pn9) = I(Pn)I(^) 

since both p„ and g are of the form (15.31) . Hence, we obtain (15.51) for / 
and g as asserted. 

Let now f,g E L2((i/i) be such that fg G h2{d'j2). Then, l{fg) is well 
defined. Let (g„) be a sequence of elements of the form (15.31) . converging 
to g in the L2((i/2) norm. Then, by the preceding argument, 

lifqn) = I(/)I(gn), Wn G N. 

In view of this equation and of the inclusion (15.41) we can write: 

||I(/^) -I(/)I((7)||l,(,,) < ||I(/c?) -I(/)I(g„)||L,(,,) + 

+ ||I(/)I(g„)-I(/)I(^?)||L,(,.) 

= Wig - qn))\\u{du)+ 
+ W)'^{<in- g)hi{du)- 
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By Cauchy-Schwartz inequality and by the isometry property of I on 
Ij2{di'), we have that 

imig-gnMuidu) < v/K(0,2vr])-||I(/(^?-g„))llL,(d.) 

= VK(0,27r])-||/(<7-g„)||L,(d^). 

Since multiphcation by / is assumed to define a bounded operator from 
L2 (ci/i) into itself, there exists a constant C > such that 

\\f{9-(ln)h2{dii)<C\\g-qn\\-L2{dji)^0 as n ^ 00. 

Using once more the Cauchy-Schwartz inequality we have: 

||I(/)I(g„ - g)\\u(du) < \mf)\\L2idu) ■ ||I(gn) - I(^)||l2H 

= ||/||L2(dM) • hn - g\\L2idji) ^0 as n ^ 00. 

The claim follows. □ 

At this stage we need a change of notation; since two (and, in the 
next section, p + 1) complex variables appear, we denote by Zi (and by 
Zi, Zp in the following section) the variables related to the Hermite 
transform, and keep z for the Z-transform variable (this notation differs 
from the one in where the Z-transform variable is denoted by C)- 

Definition 5.7. The function 

00 

(5.6) jr(z,^i) = ^z"I(Q(zi) 

n=0 

is called the generalized transfer function of the system. 

Taking the Hermite transform on both sides of (13.91) . or, equivalently, 
taking the Z transform and the Hermite transform on both sides of 
(13.71) . we obtain 

^{z, z,)=^{z,z,)^{z,z,), 

where '^{z,zi) = X^^o '^"-'-(^)('^i)' similarly for W{z,zi). The 
function Jif is analytic in a neighborhood of (0, 0) G C^. It is of 
interest to relate the properties of Jif and of the system. This is done 
in Sections 9 and 10 of the paper. We first study, in the next section, 
the case where F has a finite number of generators. 

6. The trigonometric moment problem in the poly-disc 

In [23j a solution is given to the i^-moment problem when K is a 
compact semi-algebraic set. The material is quite deep, and cannot 
be easily summarized in a short overview here. The purpose of this 
section is to serve as a guide to the reader to the topic. The starting 
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point is a semi- algebraic subset of M", defined by the positivity of k 
polynomials 

K = {x eW; pj{x) > , j = 1, . . . k} . 

Because of the application we have in mind in the next section, we will 
assume: 

Hypothesis 6.1. 

a) n is even and we set n = 2p 

b) The polynomials pj are of even degree and their highest degree ho- 
mogeneous parts have only the origin as common zero. 

One denotes by C^{K) the cone of polynomials positive on K. In [23| 
Theorem 1.4, p. 972] it is proved that, under Hypothesis 16. ![ positive 
polynomials on K belong to the additive cone 

where denotes the convex cone generated by all squares of polyno- 
mials in C[x]. The key result of [2^ is: 

Theorem 6.2. ( ]23l Lemma 3.2 p. 978]^. A functional L on M[,t] 
which is positive on C is of the form 

L{P) = [ P{x)du{x), P e R[x], 
Jk 

where dv is a positive measure on K . 

This result gives the solution to the moment problem on K: Let ai^m 
with m G Nq be complex numbers. Then there exists a positive 
measure on K such that 

/ z^z*'^dv{x) = ae^rn 
Jk 

if and only if the following conditions hold: 

L{\P{z,z*)\^)>0, VPgC[x], 
(6.1) L{p,{xMz)\')>0, \fseC[z], j = l,...p, 

L{p,{x)\P{z,z*)\^)>0, J=p+1,...K. 

7. The case of a finite number of generators 

We now assume that the Abelian group F has a finite number, say p, of 
generators, which we will denote by 71, . . . , 7p. We assume that there 
are independent in the sense that if 

7r ° ■ ■ ■ ° = ^ 
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for some integers ni, . . . G Z, then rii = ■ ■ ■ = Hp = 0. In particular, 
each generator is of the form fl2.ip . 

with 9 fixed, and where the set generates a free discrete sub- 

group of the multiphcative group of positive real numbers. We use in 
a free way the multi-index notation. 

Theorem 7.1. There is a positive measure dv on the distinguished 
boundary of the poly- disc such that 

[a{^^^)---ai^;^)di2ia)= [ e'^^'^ ■ ■ ■ e'^-'-duiO,, . . . ,9^). 
Jf Jjp 

To prove this theorem we specialize the results of the preceding section 
to the case of the poly-disc D^. It is a compact algebraic set, with 
n = 2p, m = 2p and polynomials 

Piix) = l-\z,\^,...,Ppix) = l-\zpf, 

and 

Pp+lix) = \Zi\'^ - 1, . . ■,P2p{x) = \Zp\^ - 1. 

Proof of Theorem I7.lt We define a linear form on polynomials in 
the variables 2:1, ... , Zp, z^, . . . , z* by 

Liz'^z*^) = jMli)r-^' ■ ■ ■ {<y{lp)r-^^dji{a). 

Let p be a polynomial in the variables zi, . . . , Zp^ . . . , z*. We write 
for short 

p{z,z*) =p{zi,...,Zp,zl,...,Zp). 
Let p{z, z*) = Y1 Ca,pz'^z*^. Then 
L(p(z,z*)) = 

and therefore we have 

(7.1) L{p{z, z*)) = j^p{a{^i), a(72), . . . , a(7i)*, ^(72)*, • • •) 

Since |p(z, 2;*) p is still a polynomial in z and z* , the following conditions 
hold: 

L[{l-z,z*)p{z,z*))=Q, j = l,2,...p, 
L(b(^,z*)P) >0. 
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□ 

Remark 7.2. The fact that the characters are of modulus 1 allows to 
prove (17. ip and (17.21) . It does not seem possible to relate our problems 
with another moment problem when p > 1 (for instance on the ball of 
CP). 

Definition 7.3. The Hermite transform of the element 

f{a) = Y,K<r) 

a 

%S 

l{f){z) = Y^h^z\ 

a 

Theorem 7.4. Let f G L2((i/i) be such that the operator of multipli- 
cation by f defines a bounded operator from L2 ((i/i) into itself. Then, 
for every g G L2 ((i/i).' 

(7.3) l{fg)=l{f)I{g). 

The proof is the same as for p = 1 (see the proof of Theorem 15. 6p . As 
in Definition 15.71 the function of p + 1 variables 

00 

J^{z,Zi, . ..,Zp) = ^Z"-I{hn)izi, ...,Zp) 
n=0 

is called the generalized transfer function of the system. 

8. BIBO STABILITY 

The system [TT^ will be called bounded input bounded output (BIBO) 
if there is an M > such that for every {un(7)} such that 

(8.1) sup ||M„(-)l|£2(r) < 00 

neNo 

the output is such that {yn{l)}'yer ^ hO^), n = 0, 1, . . ., and it holds 
that 

(8.2) sup ||y„(-)l|£2(r) < M sup \\un{-)\\i,(^r)- 

neNo neNo 

The following theorem gives a characterization of BIBO systems. The 
proof follows the proof of [U Theorem 3.2]. We note the following dif- 
ference between the two theorems: in [3] the multiplication operators, 
that is the counterparts of the operators Aih„ defined here using the 
Wick product, are automatically bounded. As explained there, this is 
due to Vage's inequality (see [HI Proposition 3.3.2 p. 118] and (3.1) 
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in [1] , and Section [H] below) . Here we do not have an analogue of this 
inequahty. 

Theorem 8.1. The system f ll.2p is bounded input bounded output if 
and only if the following two conditions hold: 

(a) The multiplication operators (13.51) 

Mh^- u^hn-ku, n = 0, 1,... 
are bounded from l2iX) into itself. 

(b) For all f (■) G with ||^^(-)ll^2(r) = ^ it holds that 



^•3) J2\\^lM\ur) < M. 



n=0 



Proof: That the condition (18. 3p is sufficient is readily seen. Indeed, 
take V G hiX) with ||'y(-) Il£2(r) = 1- From (II. 2p we have: 



(8.4) {Vn , v)i^(^r) = ^ (m™ , A^L-^^)^2(r) , n = 0, 1, . . . 

m=0 

and hence 

n 

\{yn, v)e^^r)\ < hm(-)ll^2(r)||-ML-™^ll^2(r) 



m=Q 



)) f E ll-^L- 

/ \m=0 



< sup \\Urn{-)\\i^(r) I I l^\\Mh„_^v\\e^(r) 

\m=0,...n 

< M sup \\um{-)\\i2(r)- 

meNo 

We obtain ([H3]) by taking v = Z/n/lbnllfeCr) when ?/„ 7^ 0. 

We now show that (18.30 is necessary. We assume that the system is 
bounded input and bounded output. We first note that the multipli- 
cation operators A4h„ are necessarily bounded. Indeed, assume that 
(18. 2p is in force and take Mq = m G ^2(r) and = for n > 0. Then, 

yn = hn-ku = Mh^u), n = 0, 1, . . . , 

and it follows from (E2D that \\MhA < M forn = 0, 1, . . . 

Let us now consider an input sequence («„) which satisfies (13. ip . For 
a given n and v choose 

Mm = if Ml M = 0, 

fl'n. — m. ' 
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and 

Ml V 
Um = -rr-r-. n Otherwise. 

We obtain from (El) and (ESD that 



m=0 

from which we get (18.31) . □ 

We now make a number of remarks: first, condition (18.21) is imphed by 
the stronger, but easier to deal with, condition 

oo 

(8.5) 5^||-M,J|<M. 

n=0 

When r is the trivial subgroup of SU{1, 1), conditions (18.31) or (18.51) 
reduce to the classical condition 

oo 

E \hn\ < OO. 
n=0 

Finally, other versions of this theorem could be given, with non causal 
systems with respect to the variable n (as in |1]), or with scale-causal 
signals. We state the last one. The proof is the same as the proof of 
Theorem 18. 1[ 

Theorem 8.2. The system (11.21) is scale-causal and hounded input 
hounded output if and only if the following two conditions hold: 
(a) The multiplication operators (13.51) 

Mh„- u^hn-ku, n = 0, 1,... 

are hounded from i2{^+) into itself. 

(6) For all v{-) G C-2iX) with \\v{-)\\(,2(v+) = I it holds that 

oo 

Ell-^L(^)ll^2(r+)<M. 

n=0 

9. DiSSIPATIVE SYSTEMS 

We will call the system (11.21) dissipative if for every input sequence (u„) 
such that 



J]ll^"(')ll^^2(r) 

n=0 



< OO 
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it holds that 



oo 



(9.1) Ell2/«(-)ll'2(r) <Ell^«("" 



^2(r)- 

n=0 n=0 



Theorem 9.1. The system is dissipative if and only if the L(£2(r))- 
valued function 



5(^) = ^^"A^,„ 



n=0 

is analytic and contractive in the open unit disc. 
Proof: Equations (19.11) expresses that the block Toeplitz operator 

(Mh, 
Mh, Mh, ■■ 

V ' ' / 

is a contraction from £2(^2(1")) into itself, and this is equivalent to the 
asserted condition on S. □ 

We consider the case of scale-causal signals (see Definition 12. 4p . 

Definition 9.2. The system (11.21) will he called scale-causal dissipative 
if the following conditions hold: 

(1) The operators M^^ are bounded from i2(^+) i'^to itself. 

(2) Condition (19. ip holds, with i2(T) replaced by£2(T^). 

Recall that we have denoted by H2((iz/) the closure in L2((iz/) of the 
powers z", where all the components of a are greater or equal to 0. 
Taking the Fourier and Hermite transforms we have: 

Theorem 9.3. The system is scale- causal dissipative if and only the 
function 

00 

(9.2) ^(^,^i,...,2p) = ^z"I(Q(zi,...,Zp) 

n=0 

is contractive from 'H.2{^) H2((iz/) into itself. Furthermore, if the 
space H2(c?i/) is a reproducing kernel Hilbert space, say with reproducing 
kernel K{zi, . . . , wi, . . .), condition (19. 2p is equivalent to the positivity 
of the kernel 

(9.3) K{zi,. . . ,wi,. . .) 

I — zw* 
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Proof: Since the operators Mh„ are assumed bounded, we have 

^n-m^Mm e £2(r+), m = 0, . . . , Tl, 

for all entries Um G £2(1"+). Thus 



and we may apply Theorem 15.61 We can write: 



/ ^ "'n—'i 



-k Ur 



n=0 



'^l{hn-m)'i-{Um) 



n=0 



H2(di^) 



Thus the dissipativity is translated into the contractivity of the block 
Toeplitz operator 

/M^ --A 
Mr Mr ■■ 

hi ho 



\ I 

from £2(H2((iz/)) into itself, and hence the claim on . To prove the 
second claim, we remark that H2(D)®H2((iz/) is the reproducing kernel 
Hilbert space with reproducing kernel 



1 — zw' 



-K{zi,--- ,wi,---). 



This comes from the fact that the reproducing kernel of a tensor prod- 
uct of reproducing kernel Hilbert spaces is the product of the repro- 
ducing kernels; see [9j , ^24j . Condition (19.31) follows then from the well- 
known characterization of bounded multipliers in reproducing kernel 
Hilbert spaces; see for instance [1], and the references therein. □ 



10. £i-i2 BOUNDED SYSTEMS 

The system f ll.2p will be called £1-^2 bounded if there is a M > such 
that for all inputs («„) satisfying 

00 

^\M-)h2{r) < 00, 

n=0 

we have 

o \ 00 

J]||2/n(-)lll(r) <MY,\M-)Uir). 

.n=0 J n=0 
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Taking the Fourier transform, this condition can be rewritten as: 

(oo \ -'-/^ OO 

n=0 J n=0 

The system (II ■2p will be called scale-causal £1-^2 bounded if it is more- 
over scale-causal, that is, if the operators Mh„ are bounded from i2(X+) 
into itself. Condition fllO.ll) then becomes: 

(00 \ 1/2 00 

n=0 / ?i=0 

from which we obtain, much in the same way as in [3], the following 
result. For completeness we present a proof. 

Theorem 10.1. A necessary and sufficient condition for the system 
(11.21) to he scalar-causal and Ix-l^ hounded is that the function 

00 

(10.3) H{z, o) = Y 2"/^(a) G H2(D) ® H2(rf/2), 

n.=0 

or, equivalently, that the transfer function 

00 

(10.4) ^{z, zi) = J2 e H2(D) ® H2{du). 

n=0 

Proof: To see that condition (110.31) is necessary, it suffices to take the 
sequence 

'1 if n = 
if n^O. 



u 

Then 



yn = hn = 0, 1, . . . , 
and condition (110.21) implies that H{z,a) G H2(D) (g) B.2{dfl). Con- 
versely, assume that the function H{z,a) G 112(115) ® H2(c?/i). From 
the expression (I3.10p . and using the Cauchy-Schwarz inequality on 

00 

Y{z, a) = H{z, a)U{z, a) = H{z, a) (z"^) , 

n=0 

we have 

00 

\\y{z,(7)\\ii2{o)<8U2{dfi) < \\H{z,a)z'''u^\\u2m<g)-H2{dfi) 

n=0 

00 

— ^ \\z'^Un\\H2{o)(8n2{d-p)\\H {z, a)\\\u2{o)(g)n2idfi)- 

n=0 
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But we have that 



\\Z Mn||H2(D)®H2(dM) — I 1 1 H2 (d/i) 5 

and so we obtain (110. 2p with M = \\H{z,a)\\\u2{o)(g)H2{d-p)- The equiv- 
alence with condition fll0.4p follows by taking the Hermite transform. 

□ 

When H.2{di') ^ h2{di') (recall that F is finitely generated), (110.41) can 
be translated into reproducing kernel conditions. In particular, in the 
cyclic case, we have: 

Theorem 10.2. Assume that H2(c?z/) ^ L2(c?z/), and let 
A{zi)A{w,)* - B{z,)B{wi)* 

1 — ZiWl 

be the reproducing kernel ofH.2{di'). The system (11.21) is scale-causal 
and ii —£2 bounded if and only if there is a M > such that the kernel 

A{zi)A{wiy -B{zi)B{wi) 



(1 - ZW*){1 - ZiWl) 

is positive in the bi-disc. 



MJ^{z,zi)je{w,wi] 



As in the case of equation (19.31) . this comes from the characterization 
of the reproducing kernel of a tensor product of reproducing kernel 
Hilbert spaces. 

11. The white noise space setting and a table 

Another kind of double convolution system, with a setting quite similar 
to the setting presented here, has been developed in [1], and rely on 
Hida's theory of the white noise space (see [19], [18], [I?] for the latter). 
We now review the main features of Hida' theory and of [3]. The 
starting point in Hida's theory is the function 



=2 II? 

"1,0 



K{si - S2) =e 

which is positive in the sense of reproducing kernels for Si,S2 in the 
Schwartz space S of real valued rapidly vanishing functions. By the 
Bochner-Minlos theorem there exists a probability measure P on the 
dual space S' of real valued tempered distributions such that 



"l2(M) 



Js' 



where we have denoted by (w, s) the duality between S and S' . The 
white noise space is defined to be the real Hilbert space L2(i5', J?-", P), 
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where T denotes the underlying Borehan sigma-algebra. 

Among all orthonormal basis of L2(iS', JF, P), there is one which plays 
a special role; it is constructed from the Hermite functions, and is 
indexed by the set £ of infinite sequences (ai,a2, • • •) indexed by N, 
and with values in No, and for which = for all j at the exception 
of at most a finite number of j. See [T9l Definition 2.2.1 p. 19]. We 
will denote by Ha (with a & i) the elements of this basis. An element 

(11.1) F = Y,UHa, faeR, 

belongs to the white noise space if 

< oo. 

The Wick product is defined by 

The white noise space is not stable under the Wick product, and there 
is the need to introduce a nuclear space, called the Kondratiev space, 
within which the Wick product is stable. The Kondratiev space is the 
projective limit of the real Hilbert spaces Tik of formal sums of the 
form flll.ip for which 

J]/^(2N)-''"<oo 

a 

for some g G N, where we use the notation 

(2N)" 2"^ X X 6"^ ■ ■ ■ 
One can also consider the complexified versions of these spaces. 

We also recall Vage's inequality (see [191 Proposition 3.3.2 p. 118]): 
Fix some integer / > 0, and let > / -|- 1. Consider h eHi and u G Hk- 
Then, h()u G Hk and 

\\h<)u\\k < A{k - l)\\h\\i\\u\\k, 

where 

A{k-l) = ^(2N)('~'=)" 

is a finite number. 

We can now introduce the systems considered in A system will be 
characterized by a sequence {hn)nez of elements in Tti for some / G N, 
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and a signal will be a sequence of elements in one of the spaces Tik-, 
with k > I Input-output relations are expressions of the form 

meZ 

Note that in view of Vage's inequality the output sequence consists 
also of elements of Tik. Furthermore, decomposing this equation along 
the basis Ha we obtain the double convolution system 

ya{n) = XI ^"-/^(^ ~ m)ufj{m), n e Z. 

meZ I3<a 

The map I which to associates the polynomial z"' is called the 
Hermite transform. It is such that 

I(fOg) = l{mg), V/,^e5_i. 

Note that under the Hermite transform the white noise space is mapped 
onto the reproducing kernel Hilbert space with reproducing kernel 
^{z,w)e2^ that is, onto the Fock space. 
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We now give the table presenting the parallels between the white noise 
space case (as applied in the paper [1]), and the present multi-scale case. 
The reader might want to look at a similar table in [3], where the analo- 
gies between the white noise space case and the hyper-holomorphic case 
are presented. 



The setting 


Stochastic case 


Multi-scale case 


Underlying 
space 


The white noise space 


^2(r) 


Hermite 
transform 


ZJ \ ~« 
i-{J^a) — Z 


I(a(7)°) = 
(F: finitely generated) 


Underlying 
reproducing 
kernel Hilbert 
space 


The Fock space 


The space H2((ii^) 


Key tool used 


Minlos theorem 
(to build the white noise 
space) 


Moment problem on the 
poly-disc (to build the 
Hermite transform) 


The product 


Wick product 


Convolution 
with respect to F. 


Double convo- 
lution 


Vain) = XlmeZ X]/3<a 

ha-i3{n - m)ui3{m) 





Remark 11.1. The pointwise product in L2 (ci/i) is a convolution in 
£2(F). Strictly speaking, it would be better to define the Hermite 
transform as the composition of the Fourier transform and of the map 
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